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SECTION  I 
INTRODUCTION 


The  transonic  compressor  has  played  a  significant  role  in 
the  development  of  modern  aircraft  propulsion  systems.  It  is  a 
major  component  in  virtually  all  modern  high-speed  aircraft 
engines  currently  in  use.  Much  time  and  effort  has  been  directed 
toward  research  into  the  fluid  dynamics  of  transonic  compressors 
to  gain  a  better  understanding  of  the  complex  flow  phenomena  which 
are  produced  by  the  rotating  machinery. 

However,  analysis  of  unsteady  flow  in  turbomachinery 
continues  to  be  one  of  the  most  difficult  and  challenging  problems 
in  internal  flows.  In  a  supersonic  or  transonic  compressor  the 
relative  Mach  number  of  the  flow  entering  the  blade  passage  is 
supersonic  while  its  axial  component  remains  subsonic.  As  a 
result  of  the  supersonic  relative  flow  entering  the  blades,  a 
complex  system  of  shock  and  expansion  waves  is  present  in  the 
blade  passage.  The  presence  of  these  waves  along  with  their 
interaction  with  the  viscous  flow  in  the  boundary  layer  results  in 
performance  degradation.  In  addition,  these  features  are 
influenced  by  the  unsteady  flow  due  to  the  wakes  from  upstream 
stators  and  an  inlet  guide  vane.  To  achieve  a  better 
understanding  of  the  interaction  between  the  unsteady  and  viscous- 
inviscid  effects  in  the  compressor,  the  use  of  the  time-dependent 
compressible  Navier-Stokes  equations  is  necessary. 

In  recent  years  significant  progress  has  been  achieved  in 
the  numerical  simulation  of  unsteady  flow  in  turbomachinery. 
Recently  Mitchell  and  Hodson  have  extended  Denton's  method J  to 
investigate  the  unsteady  inviscid  interaction  of  upstream  wakes  on 
a  moving  blade  row.  To  simulate  the  unsteady  flow  entering  the 
computational  domain  Hodson  formulated  inflow  boundary  conditions 
in  which  he  specified  the  stagnation  temperature  (assumed  constant 
in  the  wake  relative  reference  frame),  the  relative  flow  angle 
(a),  and  the  mass  flow  per  unit  area  (pu).  The  wake  is  accounted 


for  by  moving  a  velocity  profile  containing  the  wake  defect  across 
the  inflow  boundary  at  the  relative  blade  speed. 


Janssens  and  Hirsch  have  utilized  finite  element  methods  to 

treat  boundary  conditions  in  developing  a  procedure  for  computing 

viscous-inviscid  interactions  in  cascades,  accounting  for  large 

separated  boundary  layer  regions  and  correcting  for  Coriolis  and 
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centrifugal  effects  on  turbulence.  Calvert  has  also  developed  a 
viscous-inviscid  interaction  method  based  on  Denton's  time 
marching  technique  to  investigate  the  steady  flow  in  several 
transonic  compressor  cascades.  His  results  show  the  presence  of 
strong  shock  waves  and  regions  of  separated  flow  for  certain 
configurations . 

There  has  also  been  progress  toward  the  solution  of  the 

time-dependent  viscous  flow  in  turbomachinery  components.  Patched 

and  overlaid  grids  have  been  developed  for  use  in  conjunction  with 

the  thin-layer  Navier-Stokes  equations  to  simulate  very  low-speed, 

unsteady  flow  in  an  axial  turbine  configuration.^  The  current 

contract  has  been  directed  toward  the  analysis  of  time-dependent 

7 

flow  in  a  supersonic  compressor  rotor.  This  analysis  accounts 
for  the  unsteady  flow  produced  by  wakes  from  the  upstream  stator 
and  inlet  guide  vane  by  specifying  the  stagnation  pressure  profile 
at  the  inflow  boundary.  This  stagnation  pressure  profile  which  is 
obtained  from  experimental  data,  is  used  to  develop  boundary 
conditions  based  on  characteristic  variables.  With  this 
formulation  the  unsteady  behavior  occurs  in  the  form  of  an  entropy 
wave  where  entropy  is  one  of  the  characteristic  variables.  This 
formulation  also  accounts  for  the  velocity  defect  in  the  wakes 
from  the  upstream  stator  and  inlet  guide  vane.  The  downstream 
boundary  condition  is  based  on  the  concept  of  a  boundary  region 
which  incorporates  a  throat  or  choke  point  as  a  means  of 
controlling  the  back  pressure  at  the  trailing  edge  of  the  blades. 
This  approach  permits  the  use  of  the  "non-ref lective"  supersonic 
outflow  boundary  conditions  at  the  outflow  plane  while  allowing 
the  unsteady  fluctuations  to  occur  more  naturally  in  the  subsonic 
flow  at  the  trailing  edge  of  the  blades. 


2 


An  additional  objective  of  this  effort  is  to  obtain  improved 
resolution  of  the  viscous-inviscid  interaction  and  in  particular 
to  more  accurately  model  the  unsteady  behavior  on  the  suction 
surface  of  the  blades.  To  accomplish  this,  the  Baldwin-Lomax 
algebraic  eddy-viscosity  turbulence  model  is  incorporated  into  the 
computational  procedure  in  regions  near  the  blade  surfaces.  The 
implementation  of  the  turbulence  model  also  required  that  grid 
refinement  be  performed  placing  more  grid  points  in  the  boundary 
layer  regions. 

The  computations  for  the  present  work  are  performed  for  a 
stream  surface  at  approximately  90  percent  of  the  blade  span  of  a 
high-performance  transonic  compressor  rotor  being  developed  in  the 
AFWAL  Aero  Propulsion  Laboratory  at  Wright-Patterson  Air  Force 
Base.  The  blade  geometry  for  this  rotor  is  shown  in  Figure  1. 

The  computation  accounts  for  a  streamtube  contraction  of  about  38 
percent  through  the  blade  passage.  (The  current  results  are 
compared  with  those  from  previous  computations  and  qualitatively 
with  experimental  observations.) 

To  more  accurately  simulate  the  periodic  behavior  of  the 
flow  entering  the  rotor,  the  computational  domain  has  been 
expanded  to  compute  the  flow  through  four  adjacent  blades  (three 
blade  passages).  This  configuration  was  selected  because  it 
closely  approximates  the  correct  ratio  of  rotor  blades  to  upstream 
stator  blades.  In  this  particular  compressor  the  rotor  contains 
74  blades  and  the  upstream  stator  contains  54  blades.  The 
computational  arrangement  corresponds  to  the  ratio  of  four  rotor 
blades  to  three  upstream  stator  blades.  Simulation  of  the  rotor 
exactly  would  require  computing  the  flow  through  at  least  half  of 
the  rotor  (37)  blades  for  which  adequate  computer  storage  is  not 
available.  This  configuration  permits  a  more  thorough  examination 
of  mass  flow  fluctuations  which  is  the  primary  objectives  of  this 
program . 


SECTION  II 
ANALYSIS 


I 


GOVERNING  EQUATIONS 


Since  the  flow  in  a  supersonic  compressor  rotor  contains 
complex  viscous-inviscid  interactions  including  shock  waves, 
boundary  layers,  flow  separation  and  unsteady  phenomena,  a 
numerical  solution  of  the  Navier-Stokes  equations  is  required  to 
accurately  simulate  the  flow  behavior.  Thus,  the  complete  time- 


dependent  compressible  Navier-Stokes  equations  are  solved  using  an 


adaptation  of  a  computer  code  developed  by  Shang  .  Variations 


of  this  code  have  been  applied  successfully  for  internal  flows  in 


ducts  ,  thus  it  seems  appropriate  for  the  present  application. 

The  governing  equations  are  written  in  two-dimensional  rectangular 


Cartesian  coordinates  as 
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where  the  dependent  variables  are  given  by 


The  flux  vectors  are  the  Cartesian  components  of  the  continuity, 
momentum,  and  energy  equations. 
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The  components  of  the  stress  tensor  are  given  by 
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Expressed  in  this  form  minor  deviations  in  flow  area  due  to 
contraction  in  the  third  dimension  may  be  approximated. 
Sutherland's  viscosity  equation,  the  equation  of  state,  and  the 
Prandtl  number  are  specified  to  close  the  system  of  equations. 


B.  GRID  SYSTEM 

To  optimize  numerical  resolution  and  simplify  the  boundary 
conditions  the  governing  equations  are  transformed  to  a  body 
oriented  coordinate  system.  Equation  (1)  is  expressed  in  the 
transformed  coordinates  where  the  transformation  is  given  by 

C  -  &(x,y) 
n  -  n(x,y) 
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where  the  transformed  coordinate  n  »  0  line  represents  the  upper 
surface  of  the  lower  blade  and  the  n  =  JL  line  is  the  last 
row  above  the  upper  surface  of  the  upper  blade.  In  the  streamwise 
direction  the  transformed  coordinate  c  =  0  line  is  the  inflow 
boundary  and  the  £  =  KL  line  is  the  outflow  boundary.  The  two- 
dimensional  grid  is  established  between  two  adjacent  rotor  blades 
and  is  extended  upstream  to  a  point  near  the  trailing  edge  of  the 
upstream  stator  and  downstream  to  a  point  well  beyond  the  leading 
edge  of  the  trailing  stator.  However,  the  downstream  stator  is 
not  actually  included.  The  grid  boundaries  extending  upstream  of 
the  leading  edge  are  tangent  to  the  upper  surface  while  the  grid 
boundaries  extending  downstream  of  the  trailing  edge  are  tangent 
to  the  lower  surface.  The  grid  points  are  clustered  adjacent  to 
the  boundaries  which  contain  the  blade  surfaces  in  order  to 
achieve  the  necessary  resolution  of  the  flow  field  in  the  boundary 
layer  regions.  This  clustering  is  produced  by  an  exponential 
stretching  function  with  a  controlling  exponent  k  as  described  in 
Reference  8. 

A  coarse  finite  difference  mesh  for  the  rotor  cascade  is 
shown  in  Figure  2.  The  actual  grid  consists  of  60  points  in  the 
streamwise  direction  and  46  points  in  the  transverse  or  tangential 
direction.  The  grid  has  been  expanded  to  include  four  blades.  A 
coarse  mesh  for  this  configuration  is  shown  in  Figure  3- 

C.  BOUNDARY  CONDITIONS 

One  of  the  primary  objectives  of  this  program  is  to  develop 
boundary  conditions  which  will  account  for  the  time-dependent 
effects  at  the  inflow  boundary.  To  achieve  this,  different 
formulations  of  boundary  conditions  were  first  examined  for  the 
steady  flow  situation.  The  flow  enters  the  cascade  with  a 
supersonic  relative  velocity  but  has  a  subsonic  axial  component. 
The  flow  at  the  downstream  boundary  is  subsonic.  Thus,  the 
boundary  conditions  required  for  the  problem  to  be  well  posed 
consist  of  three  at  the  inflow  boundary  and  one  at  the  outflow 
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boundary .  1 1  For  this  flow  situation,  a  complex  system  of  Mach 

waves  and  shock  waves  forms  at  the  blade  leading  edge. 

Specifically,  expansion  waves  form  on  the  suction  surface  of  the 

blade  at  the  leading  edge  and  propagate  upstream  into  the  incoming 

flow.  (This  alters  the  flow  direction  and  velocity  and  the 

resulting  flow  is  treated  as  Prandtl-Meyer  flow.)  As  described  in 

References  3.  7,  12,  and  13,  the  flow  entering  the  rotor  must  be 

periodic  from  one  blade  passage  to  the  next  and  as  a  result  of  the 

expansion  waves  propagating  upstream,  a  unique  relationship  exists 

between  the  inlet  Mach  number  and  incoming  flow  angle.  This 

unique  incidence  relationship  establishes  a  unique  mass  flow  which 

must  be  accurately  predicted  since  small  errors  in  the  mass  flow 

result  in  incorrect  angles  of  attack  and  off-design  performance 

penalties.  For  these  flow  conditions,  common  practice  has  been  to 

specify  the  stagnation  pressure,  stagnation  temperature,  and 

relative  wheel  speed  at  the  inflow  boundary  while  specifying 
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static  pressure  at  the  outflow  boundary.  *  ’  In  addition  to 

these  boundary  conditions,  periodic  boundary  conditions  are 
imposed  on  the  boundaries  extending  upstream  and  downstream  of  the 
blades.  This  is  necessary  to  satisfy  the  unique  incidence 
condition  upstream  and  to  achieve  the  desired  wake  behavior 
downstream  of  the  blade  surfaces.  Finally,  the  blade  surface 
boundary  conditions  consist  of  no  slip,  and  a  prescribed  wall 
temperature.  The  compatibility  condition  is  satisfied  by  setting 
the  normal  pressure  gradient  equal  to  zero  at  the  wall.  These 
conditions  are  shown  in  Figure  4. 

A  formulation  of  the  inflow  and  outflow  boundary  conditions 
has  been  developed  based  upon  the  use  of  characteristic  variables 
in  accordance  with  the  procedure  described  by  McKenna,  Grahm  and 

Q 

Hankey  in  Reference  15.  As  Shang  notes,  the  total  number  of 
required  boundary  conditions  for  the  inflow  and  outflow  boundaries 
is  the  same  as  the  number  of  characteristic  variables  in  the 
system  of  governing  equations.  Ideally,  we  would  like  to  apply 
the  upstream  and  downstream  boundary  conditions  at  upstream  and 
downstream  infinity,  respectively.  However,  current  computer 
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storage  and  speed  limitations  dictate  that  the  boundaries  be 
placed  immediately  upstream  and  downstream  of  the  region  of 
interest.  Application  of  these  characteristic  boundary  conditions 
requires  that  they  be  placed  at  locations  where  the  assumption  of 
inviscid  inflow  is  reasonable.  In  general,  the  characteristic 
variables  are  given  by 


K1  =  p  - 


K2  =  v< 


K,  =  +  u 

3  pc 


K,  =  - - u 

11  P  c 


where  the  <*>  subscript  denotes  the  upstream  freestream  value.  The 
nature  of  the  problem  again  requires  three  boundary  conditions 
upstream.  At  the  upstream  boundary,  the  first  three 
characteristic  variables,  which  represent  vorticity,  entropy,  and 
the  downstream  running  Riemann  invariant,  are  prescribed  using 
freestream  values.  These  variables  are  then  used  in  linear 
combinations  with  the  primitive  variables  to  write  the  boundary 
conditions  in  conservative  variables. 

For  steady  state  inflow,  the  upstream  conditions  are 
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Implementation  of  these  boundary  conditions  is  intended  to 
minimize  or  eliminate  the  influence  of  reflected  waves  which 
result  from  specifying  pressure  at  inflow  or  outflow  boundaries. 
Generally,  formulations  involving  characteristic  variables  are 
developed  for  inviscid  flow  and  have  been  used  in  inviscid  flow 
calculations  for  turbomachinery . 1 1 ’ 1 ^  Since  the  present 
investigation  includes  viscous  effects,  alternative  boundary 
conditions  are  required  at  the  downstream  or  outflow  boundary,  the 
boundary  layer  and/or  blade  wake  regions. 

In  addressing  the  unsteady  flow  we  must  remember  that  a 
complex  system  of  shock  and  expansion  waves  is  established  at  the 
leading  edge  of  each  blade.  This  situation  is  complicated  further 
when  considering  the  effects  of  the  wakes  from  the  upstream  stator 
and  inlet  guide  vanes.  These  wakes  result  in  a  velocity  defect  in 
the  incoming  flow  which  changes  not  only  the  magnitude  but  also 
the  direction  of  the  incoming  velocity  profile.  Consequently,  the 
strength  and  orientation  of  the  wave  system  in  each  blade  passage 
varies  in  a  periodic  fashion  associated  with  the  blade  passage 


13 


frequency  of  the  upstream  stator  and  inlet  guide  vanes.  The 
influence  of  the  time  varying  inflow  on  the  wave  pattern  at  the 
inlet  of  the  blade  passage  is  of  particular  importance  in 
compressor  performance  since  the  inlet  wave  pattern  governs  the 
mass  flow  through  the  compressor.  We  suspect  that  the  unique 
incidence  related  to  the  circumferentially  averaged  unsteady  flow 
is  not  necessarily  the  same  as  that  related  to  an  equivalent 
steady  flow.  For  unsteady  flow  in  the  rotor,  the  characteristic 
boundary  conditions  are  applied  with  the  time-dependent  inflow 
established  by  using  the  experimentally  measured  total  pressure 
profile  to  define  entropy,  i.e.,  ,  the  first  characteristic 

variable.  The  total  pressure  in  the  stator  depicted  in  Figure  5 
is  described  by  means  of  a  Fouri er-Ser i es . 


P  (0)  -  E(A  sin  n  9  +  B  cos  n  0) 
os  n  n 


(23) 


where  0  =  ^  -  ot ;  <u  =  rotation  rate,  R  =  rotor  blade  radius. 
K 


Additionally,  the  flow  angle,  a,  exiting  the  stator  is 

assumed  invariant.  Therefore,  a  relationship  between  P  ,  u,  p 

o  s 

and  p  may  be  obtained  through  the  definition  of  total  pressure. 
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From  this  equation  density  may  be  obtained 
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and  the  velocity  defect  is  accounted  for  by 
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The  representation  of  the  characteristic  variables  used  to. 
complete  the  description  of  the  upstream  boundary  conditions  is 


P0  (  RTq) 


K1  ( 0)  =  =  — yry  ;  (entropy  profile) 
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Tnis  formulation  differs  from  that  used  in  previous  studies  in 
that  the  stator  values  of  whirl  velocity,  static  temperature  and 
total  temperature  all  vary  with  time.  Although  we  believe  that 
these  are  rational  boundary  conditions,  time  accurate  experimental 
measurements  of  these  quantities  are  required  to  confirm  the  use 
of  these  boundary  conditions. 


Downstream 

Great  difficulty  has  been  experienced  by  the  computational 
fluid  dynamics  (CFD)  community  in  computing  flows  in  short 
channels  with  subsonic  outflow.  This  is  especially  true  for 
unsteady  flows  in  which  it  is  intuitively  unnatural  to  specify  any 
flow  quantity  as  being  fixed  at  the  outflow  boundary  unless  that 
boundary  is  placed  at  downstream  infinity.  The  difficulty  was 
overcome  in  this  study  by  inserting  a  convergent-divergent  "second 
throat"  at  the  exit.  In  this  manner  supersonic  outflow  may  be 
achieved.  This  is  known  to  produce  no  numerical  difficulty.  To 
obtain  a  desired  downstream  pressure  the  "second  throat"  area  may 
be  adjusted  as  follows. 
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Recall  from  shock-relationships 


(31  ) 


With  the  implementation  of  the  "second  throat"  the  outflow 
boundary  conditions  are  simply 


where  U  is  given  in  Equation  (2). 

While  this  approach  is  still  not  a  true  simulation  of  the 
downstream  boundary,  i.e.,  stator  blade  row,  it  provides  a  more 
accurate  representation  of  the  physical  situation  than  imposing  a 
constant  variable,  such  as  static  pressure  at  the  outflow 
boundary . 

Lateral  Boundaries 

Periodic  boundary  conditions  are  imposed  upstream  and 
downstream  of  the  leading  and  trailing  edges  of  the  outermost 
blades.  At  the  blade  surfaces  no  slip,  and  a  specified  wall 
temperature  are  imposed.  These  conditions  along  with  the  inflow 
and  outflow  boundary  conditions  are  shown  in  Figure  4. 

The  evolution  of  these  boundary  conditions  is  summarized  in 
Table  1 .  The  formulation  presented  here  appears  to  give  a 
realistic  representation  of  the  time  varying  flow  effects  in  a 
transonic  compressor  rotor.  However,  it  must  be  reemphasized  that 
time  accurate  experimental  measurements  of  the  proper  flow 
quantities  are  essential  for  the  verification  of  these  boundary 
conditions . 


D.  INITIAL  CONDITIONS 

To  start  the  time-dependent  solution,  a  set  of  initial 
conditions  is  required  for  the  four  dependent  variables  at  all 
stations.  Arbitrarily,  the  upstream  conditions  were  uniformly 
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TABLE  1.  Boundary  Conditions 
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prescribed  everywhere  as  initial  conditions.  With  these 
conditions  numerical  difficulties  were  observed  at  the  downstream 
boundaries  where  periodicity  was  prescribed.  To  alleviate  this 
problem  the  periodic  boundary  conditions  were  replaced  with  solid 
surfaces  and  the  solution  permitted  to  evolve  until  the  flow 
became  aligned  with  the  channel.  After  approximately  one 
characteristic  time  period  the  correct  periodic  boundary 
condition  was  inserted  and  the  solution  procedure  continued  to 
convergence  without  further  difficulty.  Although  we  envision  many 
other  starting  procedures,  the  above  combination  of  initial 
conditions  was  used  successfully  in  all  cases  reported  herein. 

E.  TURBULENCE  MODEL 

The  turbulence  is  simulated  by  the  algebraic  eddy-viscosity 
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model  developed  by  Baldwin  and  Lomax  .  This  turbulence  model  is 
incorporated  into  the  computational  procedure  only  in  regions  r.ear 
the  blade  surfaces,  since  use  of  a  steady  flow  turbulence  model 
throughout  the  entire  computational  domain  would  tend  to 
unrealistically  damp  the  unsteady  flow  oscillations.  The 
procedure  used  for  the  present  calculations  assumes  the  presence 
of  two  layers.  The  Prandtl-Van  Driest  formulation  is  used  for  the 
inner  layer: 


e  . 
1 


pi 


(32) 


where 


l  =  ky[1  -  exp(-  y  +  /A  +  )  ] 


|w|  is  the  magnitude  of  the  vorticity: 


and 


y  ■  l  ^r> 

y  is  the  distance  normal  to  the  blade  surface,  A  +  =  26  is  the 
sublayer  thickness,  k  =  0.40  is  the  von  Karman  constant,  and  the 
subscript  w  denotes  values  at  the  blade  surfaces. 

The  model  switches  from  the  Van  Driest  formulation  to  the 
formulation  for  the  outer  region  at  the  smallest  value  of  y  for 
which  the  inner  and  outer  values  of  the  eddy-viscosity  are  equal 
(i.e.,  >  eQ).  The  formulation  for  the  outer  layer  is  given  by 

e0  =  pKCcp  Fmax  ymax  FKLEB  (33) 


where 


Fmax  '  ymax  1“  1 Cl  '  exP<-  yVs*,l 

y  is  the  value  of  y  at  which  F  occurs, 
max  max 

FKLEB  ^  +  5,5  (CKLEB  y/ymax)  ^ 
and  K  -  0.0168,  C  -  1.6,  C„. __  =  0.3. 

C  P  I\  Li  t  D 

The  turbulence  calculation  is  initiated  at  a  streamwise 

location  just  upstream  from  the  separation  point  observed  in  the 
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laminar  computations.  While  the  incorporation  of  this 
turbulence  model  should  give  improved  results,  it  is  still 
inadequate  for  accurate  prediction  of  the  flow  details.  To 
resolve  a  turbulent  boundary  layer  the  first  point  should  be 
within  the  laminar  sublayer.  An  examination  of  the  current 
results  indicated  that  the  first  point  from  the  blade  surface  was 
nominally  at  a  y+  of  30,  which  marginally  meets  this  criteria. 
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The  Baldwin-Lomax  turbulence  model  was  developed  for  steady 
flow  computations.  There  is  still  a  need  for  an  "unsteady 
turbulence"  model  which  models  only  the  fine  scale  and  does  not 
suppress  the  large  scale  fluctuations  which  are  paramount  to  this 
case . 
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SECTION  III 
NUMERICAL  PROCEDURE 


The  flow  field  being  investigated  has  a  Reynolds  number  of 

1.85x10^  based  on  a  blade  chord.  The  dominant  fluctuations  in 

this  case  are  due  to  the  interaction  of  the  rotor  with  the  stator 

wakes.  These  fluctuations  have  a  frequency  of  approximately  18.4K 

-5 

Hz,  giving  a  characteristic  time  of  5.4x10  seconds.  To  achieve 

the  temporal  resolution  for  the  high-frequency  oscillations, 

MacCormack's  explicit  and  unsplit  finite  difference  scheme  is 

appropriate  for  the  present  study.  This  algorithm  has  been  used 

successfully  for  a  wide  range  of  external  and  internal  flows 

involving  unsteady  phenomena . 1 0  The  computer  code  used  in 

this  program  is  a  degenerate  version  of  a  three-dimensional  code 
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developed  by  Shang  ’  for  operation  on  a  vector  processor.  This 
algorithm  requires  a  combination  of  alternating  forward  and 
backward  differences  for  the  predictor  and  corrector  sweeps.  This 
is  accomplished  by  using  a  different  indexing  procedure  for  the 
predictor  and  corrector  sweeps.  These  two  steps  are  defined  in 
general  by: 


U 


i.j 


U 


i.j  ~  Ax 


\  U 


i  +  1  ,  j 


At 

Ay 


» j  +  i 


(predictor ) 


(34) 


U 


n+1 
i .  J 


1  /2 


{  un  . 

I  i » J 


U 


n  +  1 

i.j 


At 

Ax 


E 


n  +  1  s 
i-  i .  y 


At  ,_n  +  1 
Ay  Q  Fi  ,  J 


F 


n  +  1 
i.j-1 


>} 


(corrector ) 


(35) 


The  efficiency  of  the  code  is  maximized  utilizing  the  CFL 

condition  on  allowable  time  increment  for  generalized  coordinates 

o 

as  reported  by  Shang': 


with  the  contravar iant  velocity  components  defined  as 


(36) 


U  -  Cxu  +  5yv 

(37) 

V  =  nxu  +  riy  v 

(38) 

The  range  of  CFL  numbers  used  in  the  present  investigation 
vary  from  0.15  to  as  high  as  0.85;  although  most  of  the  calcula¬ 
tions  have  been  made  using  a  CFL  number  of  0.8.  This  value  3eems 
to  provide  stable  results  in  a  reasonable  computational  time. 

Numerical  damping  is  also  used  in  the  present  analysis  in 
order  to  suppress  numerical  oscillations.  Again  the  form  of  the 

9 

damping  terms  is  reported  by  Shang  .  These  damping  or  artificial 
viscosity  terms  are  implemented  in  each  sweep  direction.  They 
are: 
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In  the  present  analysis  the  value  of  the  damping  constant  (ft)  used 
is  2. 


For  steady  inflow  conditions  the  computations  were  performed 
for  6,400  time  steps  which  corresponds  to  about  8.9  characteristic 
times  based  on  chord  length  (L)  divided  by  the  upstream  freestream 
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velocity.  This  required  about  3.5  hours  of  CYBER  845  computer 
time.  The  computations  using  the  time-dependent  inflow  conditions 
were  carried  out  to  6,000  time  steps  at  which  point  the  turbulence 
model  was  implemented  and  the  computation  was  continued  to  11,500 
time  steps  to  ensure  that  a  limit  cycle  was  realized.  This 
corresponds  to  about  16  characteristic  times  based  on  (L/U^). 


24 


SECTION  IV 

DISCUSSION  OF  RESULTS 


In  this  investigation  computations  are  initially  performed 
for  flow  through  two  adjacent  blades  in  a  transonic  compressor 
rotor,  at  a  radial  location  corresponding  to  approximately  90 
percent  of  the  span.  This  rotor  has  a  relative  inflow  angle  of  61 
degrees  from  the  axial  direction.  This  results  in  a  subsonic 
axial  velocity  component  having  a  Mach  number  of  0.66. 

The  primary  goal  of  this  program  is  to  examine  and  assess 
the  influence  of  unsteady  flow  effects  on  the  mass  flow  through 
the  compressor  rotor.  Two  significant  sources  of  unsteady  flow 
are  (1)  the  wakes  from  the  upstream  stator  and  inlet  guide  vane 
and  (2)  the  oscillatory  flow  separation  and  vortex  shedding  on  the 
suction  surface  of  the  blade.  The  wakes  from  the  upstream  blade 
rows  enter  the  rotor  blade  passage  at  a  well  defined  periodicity 
determined  by  the  blade  passage  frequency  of  the  rotor.  The 
resulting  unsteady  flow  is  classified  as  a  forced  oscillation.  In 
contrast  the  extent  of  the  flow  separation  on  the  blade  suction 
surface  varies  with  the  shedding  or  vortices  from  the  blade 
trailing  edge.  This  vortex  shedding  does  not  occur  with  the  same 
frequency  or  period  between  subsequent  vortices  hence  this 
unsteady  behavior  is  quasi-periodic  and  is  classified  as  a  self- 
sustained  oscillation. 

Both  of  these  unsteady  phenomena  affect  the  mass  flow  rate. 

In  order  to  analyze  the  influence  of  these  phenomena, 

computations  were  first  performed  for  steady  inflow  boundary 

conditions  thereby  eliminating  the  wakes  from  the  upstream  blade 

rows.  These  computations  represent  laminar  flow  through  single 

blade  passage  with  attention  focused  on  separated  flow  on  the 

blade  suction  surface.  Typical  features  of  the  flow  separation 

with  the  steady  inflow  boundary  conditions  are  shown  in  Figure  6. 

The  presence  of  separated  flow  in  such  regions  has  been  noted  by 
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other  researchers  in  previous  efforts  ’  .  However,  those 
investigations  treated  only  steady  flow. 


Since  the  primary  interest  of  the  present  program  is  on  the 
influence  of  unsteady  behavior  on  mass  flow,  the  upstream  boundary 
conditions  were  modified  to  account  for  the  unsteady  flow  due  to 
the  wakes  from  the  upstream  blade  rows.  Computations  were  carried 
out  for  several  characteristic  times  assuming  purely  laminar 
unsteady  flow.  The  mass  flow  rate  for  this  configuration  is  shown 
in  Figure  7.  This  figure  shows  the  periodic  fluctuation  of  the 
mass  flow  which  results  from  the  wakes  from  the  upstream  blade 
rows.  The  unsteady  mass  flow  rate  also  indicates  that  a  secondary 
oscillation  of  much  lower  frequency  is  present.  This  flow 
frequency  oscillation  seems  to  be  associated  with  the  unsteady 
behavior  of  the  region  of  separated  flow  on  the  blade  suction 
surfaces.  The  quasi-periodic  nature  of  the  flow  separation  i3 
related  to  the  vortex  shedding  from  the  blade.  The  physical 
interpretation  of  this  phenomenon  is  that  the  vortex  forms  on  the 
blade  suction  surface  and  grows  as  it  migrates  along  the  surface 
in  the  streamwise  direction.  This  results  in  growth  of  the 
separation  which  acts  as  a  blockage  in  the  blade  passage.  This 
blockage  then  reduces  the  mass  flow  through  the  blade  passage  in 
an  oscillatory  fashion  corresponding  to  the  vortex  shedding 
frequency.  Figure  7  also  shows  a  comparison  between  the  time 
averaged  non-steady  mass  flow  rate  and  the  circumferentially  mass- 
averaged  uniform  inflow  into  the  rotor.  It  is  significant  that 
the  time  averaged  non-steady  mass  flow  is  about  5  percent  higher 
than  the  mass  averaged  uniform  inflow. 

The  magnitude  of  the  flow  separation  is  much  larger  than 
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that  observed  in  previous  studies  ,  as  is  the  amplitude  of  the 
mass  flow  fluctuation  due  to  the  wakes  entering  the  blade  passage. 
Thus,  it  is  believed  that  the  influence  of  turbulence  effects  are 
not  adequately  taken  into  consideration  with  this  formulation. 
Therefore,  a  turbulence  model  was  incorporated  into  the 
computational  procedure  in  regions  near  the  blade  surfaces  in 
order  to  account  for  these  effects.  The  turbulence  model  used  in 
these  computations  is  the  algebraic  eddy  viscosity  model  of 


[ 

! 
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1  Baldwin  and  Lomax.  The  incorporation  of  any  turbulence  model  must 

be  accomplished  with  extreme  caution  since  they  were  developed  for 
steady  flow.  Hence,  improper  implementation  can  result  in 
suppression  of  the  unsteady  flow  features. 

The  oscillatory  mass  flow  for  the  turbulent  computation  is 
shown  in  Figure  8.  In  comparing  this  figure  with  the  laminar 
results  (Figure  7),  it  is  seen  that  the  amplitude  of  the 
oscillation  due  to  the  incoming  wakes  has  been  reduced 
substantially.  It  also  appears  that  the  low  frequency  oscillation 
associated  with  the  vortex  shedding  has  been  suppressed  as  well. 
Also  in  this  figure  the  time-averaged  non-steady  mass  flow  rate  is 
compared  with  the  mass  flow  computed  for  a  uniform  inflow 
corresponding  to  the  circumferentially  mass-averaged  rotor  inflow. 
Note  that  as  in  the  laminar  case  the  time  averaged  non-steady  mass 
flow  rate  is  higher  than  the  corresponding  thermodynamically 
equivalent  steady  flow,  although  the  difference  is  only  about 
three  percent  when  the  turbulence  model  is  included  compared  to  5 
percent  in  the  laminar  case. 

Further  insight  can  be  gained  by  examining  other  flow 
properties.  The  extent  of  the  separation  for  the  laminar 
computation  with  unsteady  inflow  boundary  conditions  is  shown  in 
velocity  vectors  in  Figure  9.  A  careful  examination  of  this  flow 
separation  reveals  the  formation  of  the  vortex  on  the  suction 
surface.  This  vortex  is  ultimately  shed  and  migrates  through  the 
blade  wake  region  in  the  direction  opposite  to  the  rotation  of  the 
blades.  Figure  10  shows  velocity  vectors  for  the  unsteady  flow 
j  using  the  Baldwin-Lomax  turbulence  model  in  the  computation. 

|  These  velocity  vectors  show  that  the  incorporation  of  the 

i 

i  turbulence  model  diminishes  the  extent  of  the  flow  separation 

i  significantly.  In  fact,  the  region  which  was  separated  in  the 

i 

I  laminar  computation  has  become  primarily  a  region  of  very  low- 

speed  flow  with  only  a  small  separation  occurring  very  near  the 
1  trailing  edge  of  the  blade  on  the  suction  surface.  Although  much 

smaller  than  in  the  laminar  case,  a  vortex  is  shed  from  the 
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’'igure  7.  Comparison  of  time-dependent  mass  flow  and  stead; 

state  mass  flow  at  blade  leading  edge  -  laminar 
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Figure  8.  Comparison  of  time-dependent  mass  flow  and  steady 
state  mass  flow  at  blade  leading  edge  -  turbulent 
calculation. 


trailing  edge  and  migrates  through  the  wake  region  in  much  the 
same  fashion  as  for  the  laminar  computation. 

The  Mach  contours  also  display  a  number  of  the  unsteady  flow 
features.  From  the  results  of  Reference  7,  the  Mach  number 
contours  show  the  oscillatory  behavior  of  the  oblique  shock  wave 
at  the  blade  leading  edge  (Figure  11).  Also  this  figure  shows 
that  the  shock  wave  impinges  on  the  suction  surface  boundary  layer 
in  the  region  where  the  onset  of  separation  seems  to  occur.  The 
oscillatory  growth  and  decay  of  the  separation  near  the  blade 
trailing  edge  is  also  seen  in  this  figure.  Comparison  with  the 
velocity  vector  plots  shows  this  behavior  is  associated  with  the 
formation  and  shedding  of  the  vortex  structure.  The  Mach  number 
contours  for  the  computations  with  turbulence  included  are  shown 
in  Figure  12.  The  time-dependent  behavior  shown  in  this  figure 
indicates  that  the  turbulence  model  reduces  the  strength  of  the 
leading-edge  shock  and  the  extent  of  flow  separation.  In 
addition,  note  that  the  oscillatory  nature  of  the  flow  separation 
is  reduced  significantly  in  amplitude  by  including  the  turbulence 
model.  The  large  fluctuations  in  the  flow  for  the  laminar  case 
are  shown  in  the  time  varying  blade  surface  pressure  distributions 
(Figure  13)-  Again  the  large  amplitude  in  pressure  fluctuation  at 
the  blade  trailing  edge  is  associated  with  the  vortex  shedding 
phenomena.  This  figure  also  reveals  a  large  spike  in  the  static 
pressure  on  both  the  pressure  and  suction  surfaces  at  the  blade 
leading  edge.  This  is  the  result  of  numerically  smearing  the 
shock  so  that  it  behaves  as  a  detached  bow  shock  in  which  the 
pressure  tends  toward  the  stagnation  value  on  the  pressure 
surface.  The  flow  then  rapidly  expands  arouna  the  leading  edge  on 
both  surfaces  and  the  pressure  recovers  toward  the  value  which 
would  be  produced  by  a  weak  attached  shock  on  the  pressure 
surface.  This  phenomena  is  smeared  over  a  much  larger  region  than 
would  be  expected  in  an  actual  physical  flow  due  to  the  relatively 
coarse  grid  in  the  vicinity  of  the  blade  leading  edge.  This 
oblique  shock  wave  then  impinges  on  the  suction  surface  of  the 
adjacent  blade  causing  boundary  layer  separation  in  a  classical 
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Figure  11.  Mach  contours  for 

laminar  computation . 


Figure  12.  Mach  contours  for 

turbulent  computation 


shock  wave--'ooundary  layer  interaction.  The  wave  is  then 
reflected  back  toward  the  pressure  surface  causing  a  subsequent 
rise  in  the  surface  pressure.  The  drop  in  static  pressure  on  the 
suction  surface  shows  the  presence  of  the  expansion  just  inside 
the  blade  passage.  Again  comparing  with  the  surface  static 
pressure  distribution  for  the  turbulent  calculation  we  observed 
that  the  magnitude  of  the  fluctuations  are  reduced  significantly 
(Figure  14).  in  particular,  the  large  fluctuation  at  the  trailing 
edge  for  the  laminar  case  is  almost  completely  eliminated. 
Referring  to  the  velocity  vectors  (Figures  9  and  10),  note  that  in 
the  turbulent  case  the  vortex  strength  at  the  trailing  edge  has 
been  reduced  to  nearly  zero  in  the  turbulent  computations.  Also 
note  that  the  sharp  spike  near  the  leading  edge  of  the  pressure 
surface  has  been  eliminated,  supporting  the  previous  observation 
that  the  shock  strength  had  been  reduced  in  the  turbulent 
calculat ions . 

Examination  of  the  stagnation  pressure  contours  for  both  the 
laminar  (Figure  15)  and  turbulent  (Figure  16)  computations  is 
revealing.  For  the  laminar  case  the  periodic  behavior  repeats  at 
the  blade  leading  edge  at  approximately  0.55  characteristic  time 
while  the  behavior  at  the  trailing  edge  is  not  repeatable  due  to 
the  oscillatory  behavior  of  the  flow  separation.  In  the  case  of 
the  turbulent  flow  the  periodic  behavior  at  the‘N  leadi  ng  edge  also 
repeats  at  about  0.55  characteristic  time.  However,  note  for  this 
case  that  this  periodicity  is  essentially  preserved  throughout  the 
blade  passage.  This  is  attributed  to  the  fact  that  the  turbulence 
model  suppresses  the  oscillatory  behavior  of  the  flow  separation 
region.  The  time  variation  of  stagnation  pressure  can  be  regarded 
as  entropy  waves  migrating  through  the  blade  passage.  Viewing 
the  stagnation  pressure  in  this  fashion  gives  a  good  view  of  the 
time  history  of  the  wakes  from  the  upstream  components. 

Computation  for  Multiple  Blade _ Passages 

Although  the  computation  of  time-dependent  flow  through  a 
single  blade  passage  provides  significant  insight  into  many 
aspects  of  the  unsteady  flow  in  the  compressor  rotor,  this 


simulation  does  not  adequately  represent  the  details  of  the 
periodic  nature  of  the  flow.  In  order  to  simulate  the  periodic 
behavior  exactly  it  would  be  necessary  to  compute  the  flow  through 
at  least  half  of  the  rotor  (i.e.,  37  blades)  since  the  ratio  of 
stator  to  rotor  blades  is  5*1:74.  Since  current  computer  storage 
is  not  sufficient  to  accommodate  a  time  accurate  computation  of 
this  magnitude  the  real  configuration  is  approximated  by  computing 
the  flow  through  three  blade  passages,  assuming  three  upstream 
stator  blades  for  every  four  rotor  blades.  The  computed  results 
for  this  configuration  much  more  closely  simulate  the  actual 
compressor  flow  than  the  computations  for  the  single  blade 
passage . 

The  computation  for  three  blade  passages  was  accomplished  by 
first  expanding  the  grid  to  this  configuration  and  then  performing 
the  calculations  for  laminar  flow  with  the  steady  inflow  boundary 
conditions.  After  a  stable  solution  was  obtained,  the  time- 
dependent  inflow  boundary  conditions  were  imposed  to  account  for 
the  wakes  from  upstream  components.  The  results  for  this 
computation  give  a  better  indication  of  how  time  dependent 
features  vary  between  blade  passages  (than  in  the  single  blade 
passage  case).  Figures  17a,  b,  c,  and  d  represent  typical 
examples  of  Mach  number,  static  pressure,  stagnation  pressure  and 
velocity  vectors  for  the  laminar  steady  in  flow  case.  It  is 
believed  that  more  accurate  simulations  are  achieved  by 
incorporating  the  Baldwin  Lomax  turbulence  model.  The  effect  of 
this  again  was  suppression  of  the  magnitude  of  the  unsteady  flow 
features . 

To  examine  the  influence  of  the  unsteady  flow  separation  on 
the  blade  suction  surface,  computations  were  performed  for  several 
characteristic  times  with  steady  inflow  boundary  conditions. 

Figure  18  shows  the  time  varying  mass  flow  at  the  blade  leading 
edge  for  this  case.  The  oscillatory  behavior  of  the  mass  flow 
rate  is  associated  with  the  quasi-periodic  fluctuations  in  the 
flow  separation  on  the  blade  suction  surface.  As  noted  previously 
this  is  related  to  the  vortex  shedding  at  the  blade  trailing  edge. 
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Careful  examination  of  this  oscillatory  flow  separation  reveals 
that  it  migrates  from  one  blade  row  to  the  next.  Although  this 
behavior  is  somewhat  reminiscent  of  rotating  stall,  it  occurs  at 
a  much  higher  frequency  than  would  be  expected  in  an  actual 
rotor.  This  can  be  attributed  to  fact  that  the  periodic  boundary 
conditions  are  imposed  over  only  four  blades.  Consequently,  the 
disturbance  can  only  migrate  through  three  blade  passages  before 
it  repeats  while  in  the  actual  rotor  it  may  travel  around  the 
full  circumference.  Since  the  flow  separation  is  strongly 
coupled  to  the  adverse  pressure  gradient  across  the  blade  row, 
the  magnitude  of  this  separation  and  the  resultant  flow  blockage 
were  reduced  by  decreasing  the  pressure  ratio  across  the  blade 
row.  This  was  accomplished  by  increasing  the  throat  area  of  the 
downstream  choke  nozzle  through  a  series  of  numerical 
experiments.  From  the  numerical  experiment  it  was  found  that 
increasing  the  downstream  nozzle  area  to  correspond  to  an  area 
ratio  of  A/A*,,  ~  0.65  eliminates  the  migration  of  the  flow 
separation  from  one  blade  passage  to  another.  The  influence  of 
varying  the  nozzle  area  ratio  on  the  mass  flow  rate  is  shown  in 
Figure  19.  The  variation  of  absolute  stagnation  pressure  ratio 
across  the  rotor  as  a  function  of  nozzle  area  is  shown  in  Figure 
20.  As  indicated  in  these  figures  the  mass  flow  rate  and  the 
stagnation  pressure  ratio  behave  somewhat  erratically  for  nozzle 
area  ratios  at  which  the  unsteady  separation  is  present.  Figure 
21  shows  the  variation  of  absolute  stagnation  pressure  with  mass 
flow  rate.  The  results  of  the  numerical  experiment  are 
summarized  in  Table  2. 

Upon  elimination  of  the  unsteady  flow  separation,  the  time- 
dependent  inflow  boundary  conditions  were  once  again  incorporated 
into  the  computation.  The  mass  flow  rate  at  the  blade  leading 
edge  is  shown  both  with  these  boundary  conditions  and  with  the 
steady  inflow  boundary  conditions  in  Figure  22.  For  the  steady 
inflow,  the  mass  flow  rate  is  0.362  slugs/ ( sec-ft )  while  the 
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Figure  20.  Computed  stagnation  pressure  ratio  as  a  function 
of  exit  nozzle  area  ratio. 
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average  mass  flow  rate  for  the  unsteady  inflow  is  0.364 
slugs/ ( sec-ft ) .  Although  the  unsteady  value  is  slightly  higher 
(less  than  1$),  the  difference  between  the  two  is  much  smaller 
than  that  observed  for  higher  pressure  ratios  across  the  blade 
row.  For  higher  pressure  ratios  the  unsteady  inflow  produced  an 
average  mass  flow  in  excess  of  3%  higher  than  that  for  steady 
inflow.  The  reduction  in  this  difference  with  pressure  ratio 
implies  that  the  oscillatory  unsteady  inflow  is  not  the  sole 
source  of  the  higher  mass  flow  rate.  It  should  also  be  noted 
that  the  overall  average  level  of  the  mass  flow  increases  as  the 
absolute  stagnation  pressure  ratio  across  the  blade  row  is 
decreased.  Additional  details  of  the  flow  for  this  configuration 
are  seen  by  examining  the  Mach  number  contours,  static  and 
stagnation  pressure  contours  and  velocity  vectors  as  they  vary 
with  time.  These  flow  properties  are  shown  in  Figures  23,  24, 

25,  and  26,  respectively. 
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SECTION  V 

CONCLUSIONS  AND  RECOMMENDATIONS 


Procedures  have  been  developed  for  treating  boundary  condi¬ 
tions  for  a  supersonic  flow  entering  a  two-dimensional  compressor 
rotor  blade  row  and  a  subsonic  flow  exiting  that  blade  row.  The 
upstream  boundary  conditions  are  formulated  using  characteristic 
variables  and  an  experimentally  measured  stagnation  pressure 
profile  to  simulate  the  pressure  and  velocity  defects  from  the 
wake  of  upstream  stators  and  inlet  guide  vanes.  This  formulation 
also  produces  the  proper  entropy  variation  at  the  inflow  boundary. 

The  downstream  boundary  conditions  are  developed  by  incorporating 

£ 

a  "second  throat"  (A2)  downstream  of  the  blade  trailing  edge  and 
imposing  supersonic  outflow  boundary  conditions  at  the  exit  of  the 
throat.  The  area  of  the  "second  throat"  can  be  adjusted  to 
relocate  the  shock  and  the  pressure  rise  across  the  blade  passage. 
The  computation  also  incorporates  streamtube  contraction  through 
the  axial  component  of  the  continuity  equation  to  allow  for  some 
three-dimensional  effects.  These  boundary  conditions  are  used  to 
obtain  solutions  to  the  complete  time-dependent  compressible 
Navi er-Stokes  equations  for  the  transonic  compressor  rotor.  A 
series  of  four  rotor  blades  with  three  stator  wakes  are  computed 
with  periodicity  imposed.  This  represents  the  periodic  behavior 
of  a  stage  with  74  rotor  blades  and  54  stator  blades. 

Both  laminar  and  turbulent  cases  are  computed.  Results  are 
presented  in  the  form  of  velocity  vectors,  Mach  and  total  pressure 
contours,  blade  surface  pressure  distribution  and  mass  flow 
history.  The  computed  results  give  a  good  representation  of  the 
time  varying  flow  properties  in  a  compressor  rotor  and  provide  a 
significant  advance  in  solving  the  time-dependent  viscous  flow 
field  in  turbomachinery. 
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Recommendations  are: 

a.  Decrease  region  of  bow  shock  influence. 

b.  Resolve  shock  reflection  region. 

c.  Improve  turbulence  model  for  unsteady  cases 

d.  Investigate  three-dimensional  effects. 

e.  Examine  rotor-stator  interaction. 
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NOMENCLATURE 


c 

e 

E,F 
i .  J 
k 

K. 

1 

L 

n 

P 

R 

t 

T 

U 


u,  v 


w 

X 
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a 

6 

Y 

6 

S  ,n 
0 
A 


y 

p 


o 

T 

0) 


XX 

xy 


speed  of  sound 
specific  energy 

vector  fluxes  in  mean  flow  equations 
indices  for  grid  point  location 
thermal  conductivity 
characteristic  variables 
characteristic  length  (blade  chord) 
time  index 
pressure 

gas  constant,  radius 
time 

temperature 
dependent  variable 
velocity  components 
wheel  speed 
axial  coordinate 
tangential  coordinate 
stator  flow  angle 
damping  factor 
ratio  of  specific  heats 
streamtube  contraction  factor 
transformed  coordinates 
azimuthal  angle 

-2/3  y,  second  viscosity  coefficient 
molecular  viscosity  coefficient 
density 

normal  stresses 
shear  stress 
vorti ci ty 
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NOMENCLATURE  (Continued) 


Subscripts 


00 


w 

JL 

r 

KL 


max 


min 


o 

s 

u 


free  stream  (annulus)  condition 

wall  condition 

maximum  value 

minimum  value 

lower  surface 

stagnation  value 

stator  coordinate  frame 

upper  surface 
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